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ABSTRACT 


The  objective  of  this  study  was  to  examine  the  relation- 
ship between  forward  ship  motion  and  the  effects  of  aero- 
dynamic and  hydrodynamic  disturbances  on  motions  in  the 
horizontal  plane.   It  includes  course  control  and  stability 
analysis  for  the  unsteered  and  steered  cases  (manual  and 
automatic)  using  several  sets  of  operating  conditions. 

Hydrodynamic  and  aerodynamic  effects  are  considered  to 
be  functions  of  hull  motion,  rudder  (steered  cases) ,  pro- 
peller and  wind  effects  ratio.   A  dimensionless  mathematical 
model  has  been  developed  and  solved  with  respect  to  ship 
axes. 


TABLE  OF  CONTENTS 

Page 

I.  INTRODUCTION 13 

II.  MATHEMATICAL  MODEL  15 

A.   COORDINATE  SYSTEMS  AND  SHIP  MOTION  ....  15 

1.   Equations  of  Motion  (on  the 

horizontal  plane)   15 

III.  CONTROL  STUDIES 27 

A.  MANUALLY  STEERED  SHIP  MOTION  31 

B.  AUTOMATICALLY  SHEERED  SHIP 37 

IV.  STABILITY  ANALYSIS  43 

V.  DISCUSSION  OF  RESULTS  AND  CONCLUSIONS   ....  54 

APPENDIX 56 

COMPUTER  OUTPUTS 103 

COMPUTER  PROGRAMS   Ill 

REFERENCES 121 

INITIAL  DISTRIBUTION  LIST   122 

FORM  DD  1473 123 


LIST  OF  TABLES 

1.  Variation  of  Yaw  for  Unsteered  Ship 33 

2.  Variation  of  Sway  for  Unsteered  Ship 34 

3.  Variation  of  Yaw  for  Manually  Steered  Ship   ...  36 

4.  Variation  of  Sway  for  Manually  Steered  Ship    .  .  36 

5.  Yaw  Gain  Constant  and  Yaw  Rate  Gain  Con- 
stant Values  for  Several  Wind  Directions 

and  Velocities 40 

6.  Variation  of  Yaw  for  Automatically  Steered 

Ship 42 

7.  Variation  of  Sway  for  Automatically  Steered 

Ship 42 

8.  Real  Roots  of  Characteristic  Equation  for 

Head  Wind 49 

9.  Real  Roots  of  Characteristic  Equation  for 

Beam  Wind 50 

10.  Real  Roots  of  Characteristic  Equation  for 

Stern  Wind 51 

11.  Real  Roots  of  Characteristic  Equation  for 

Wind  Velocity  =  7  Miles 52 

12.  Real  Roots  of  Characteristic  Equation  for 

Wind  Velocity  =  9  Miles 53 


LIST  OF  FIGURES 

Page 

1.  Relation  of  Wind  Velocity  and  Direction 

to  the  Ship  Axes 19 

2.  Positive  Direction  of  Yaw 20 

3.  Hydrodynamic  Coefficients  Versus  Rudder 

Force  Coefficients 25 

4.  Computer  Model  for  Program  II  and  III 28 

5.  Variation  of  Yaw  Versus  Time  for  Wind 

Direction  000° ,  Velocity  1  mile/hr  .  .     58 

6.  Variation  of  Sway  Versus  Time  for  Wind 

Direction  000° ,  Velocity  1  mile/hr  59 

7.  Variation  of  Yaw  Versus  Time  for  Wind 

Direction  000°,  Velocity  3  miles/hr   60 

8.  Variation  of  Sway  Versus  Time  for  Wind 

Direction  000°,  Velocity  3  miles/hr   61 

9.  Variation  of  Yaw  Versus  Time  for  Wind 

Direction  000° ,  Velocity  5  miles/hr   62 

10.  Variation  of  Sway  Versus  Time  for  Wind 

Direction  000° ,  Velocity  5  miles/hr   63 

11.  Variation  of  Yaw  Versus  Time  for  Wind 

Direction  090° ,  Velocity  3  miles/hr   64 

12.  Variation  of  Sway  Versus  Time  for  Wind 

Direction  090° ,  Velocity  3  miles/hr   65 

13.  Variation  of  Yaw  Versus  Time  for  Wind 

Direction  090°,  Velocity  5  miles/hr   66 

14.  Variation  of  Sway  Versus  Time  for  Wind 

Direction  090° ,  Velocity  5  miles/hr   67 

15.  Variation  of  Yaw  Versus  Time  for  Wind 

Direction  090°  ,  Velocity  7  miles/hr   68 

16.  Variation  of  Sway  Versus  Time  for  Wind 

Direction  090° ,  Velocity  7  miles/hr   69 

17.  Variation  of  Yaw  Versus  Time  for  Wind 

Direction  180° ,  Velocity  1  mile/hr  70 

18.  Variation  of  Sway  Versus  Time  for  Wind 

Direction  180° ,  Velocity  1  mile/hr  71 


19.  Variation  of  Yaw  Versus  Time  for  Wind 

Direction  180°,  Velocity  3  miles/hr  ......      72 

20.  Variation  of  Sway  Versus  Time  for  Wind 

Direction  180° ,  Velocity  3  miles/hr 73 

21.  Variation  of  Yaw  Versus  Time  for  Wind 

Direction  180° ,  Velocity  5  miles/hr  74 

22.  Variation  of  Sway  Versus  Time  for  Wind 

Direction  180°,  Velocity  5  miles/hr  75 

23.  Variation  of  Yaw  Versus  Time  for  Wind 

Direction  180° ,  Velocity  7  miles/hr  76 

24.  Variation  of  Sway  Versus  Time  for  Wind 

Direction  180° ,  Velocity  7  miles/hr 77 

25.  Variation  of  Yaw  Versus  Time  for  Wind 
Direction  090°,  Velocity  5  miles/hr, 

Rudder  Deflection  15°  78 

26.  Variation  of  Sway  Versus  Time  for  Wind 
Direction  090°,  Velocity  5  miles/hr, 

Rudder  Deflection  15°  79 

27.  Variation  of  Yaw  Versus  Time  for  Wind 
Direction  180°,  Velocity  3  miles/hr, 

Rudder  Deflection  15°  80 

28.  Variation  of  Sway  Versus  Time  for  Wind 
Direction  180°,  Velocity  3  miles/hr, 

Rudder  Deflection  15°  81 

29.  Variation  of  Yaw  Versus  Time  for  Wind 
Direction  180°,  Velocity  5  miles/hr, 

Rudder  Deflection  15°  82 

30.  Variation  of  Sway  Versus  Time  for  Wind 
Direction  180°,  Velocity  5  miles/hr, 

Rudder  Deflection  15°  83 

31.  Variation  of  Yaw  Versus  Time  for  Wind 
Direction  000°,  Velocity  1  mile/hr, 
Automatically  Steered  84 

32.  Variation  of  Sway  Versus  Time  for  Wind 
Direction  000°,  Velocity  1  mile/hr, 
Automatically  Steered  85 

33.  Variation  of  Yaw  Versus  Time  for  Wind 
Direction  000°,  Velocity  3  miles/hr, 
Automatically  Steered  86 


34.  Variation  of  Sway  Versus  Time  for  Wind 
Direction  000°,  Velocity  3  miles/hr, 
Automatically  Steered  87 

35.  Variation  of  Yaw  Versus  Time  for  Wind 
Direction  000°,  Velocity  5  miles/hr, 
Automatically  Steered  88 

36.  Variation  of  Yaw  Versus  Time  for  Wind 
Direction  090° ,  Velocity  3  miles/hr, 
Automatically  Steered  ....      89 

37.  Variation  of  Sway  Versus  Time  for  Wind 
Direction  090°,  Velocity  3  miles/hr, 
Automatically  Steered  90 

38.  Variation  of  Yaw  Versus  Time  for  Wind 
Direction  180°,  Velocity  1  mile/hr, 
Automatically  Steered  91 

39.  Variation  of  Sway  Versus  Time  for  Wind 
Direction  180°,  Velocity  1  mile/hr, 
Atuomatically  STeered  92 

40.  Yaw  Gain  Constant  Versus  Yaw  Rate  Gain 
Constant  for  Wind  Direction  000°, 

Velocity  1  mile/hr 93 

41.  Yaw  Gain  Constant  Versus  Yaw  Rate  Gain 
Constant  for  Wind  Direction  000°, 

Velocity  3  miles/hr  94 

42.  Yaw  Gain  Constant  Versus  Yaw  Rate  Gain 
Constant  for  Wind  Direction  000°, 

Velocity  5  miles/hr  95 

43.  Yaw  Gain  Constant  Versus  Yaw  Rate  Gain 
Constant  for  Wind  Direction  090°, 

Velocity  3  miles/hr  96 

44.  Yaw  Gain  Constant  Versus  Yaw  Rate  Gain 
Constant  for  Wind  Direction  180°, 

Velocity  1  mile/hr 97 

45.  Real  a.  Versus  Wind  Velocity  for  Head  Wind   .  .     98 

46.  Real  a..  Versus  Wind  Velocity  for  Beam  Wind   .  .     99 

47.  Real  a.  Versus  Wind  Velocity  for  Stern  Wind  .  '.     100 

48.  Real  a.  Versus  Wind  Direction  for  7  miles/hr 

Wind  Velocity 101 

49.  Real  a.  Versus  Wind  Direction  for  9  miles/hr 

Wind  velocity 102 

7 


PHYSICAL  QUANTITIES 


a.  Coefficients  of  perturbation  equation  in  yaw 

b.  Coefficients  of  perturbation  equation  in  sway 

c.  Coefficients  of  perturbation  equation  in  surge 
F  Force  on  the  body 

J    Moment  of  inertia  referred  to  z-axis 

z 

I  Ship  length 

M    Moment  about  origin 

m    Mass  of  ship 

N    Hydrodynamic  and  aerodynamic  yaw  moment 

N    Aerodynamic  yaw  moment 

N«    Derivative  of  hydrodynamic  yaw  moment  with  respect  to 

yaw  acceleration 
N    Derivative  of  hydrodynamic  yaw  moment  with  respect  to 

sideslip  velocity 
N    Derivative  of  hydrodynamic  yaw  moment  with  respect  to 

rudder  angle 
R    Vector  distance  of  a  point  from  the  origin,  and  the 

components  along  the  X,  Y,  and  Z  axis   R  =  ix  +  jy  +  kz 

r    Yaw  rate 

-*■■■?■    -t  ■* 

U    Velocity  of  the  origin  U  =  iu  +  jv  +   kw 


u    Component  of  equilibrium  ship  speed  in  x-axis  direction 

u    Component  of  ship  velocity  relative  to  air  in  x-axis 

direction 

U    Ship  velocity  relative  to  air 
a 

U    Wind  velocity  relative  to  earth 


v    Component  of  ship  speed  in  y-axis  direction 

v    Component  of  ship  velocity  relative  to  air  in  y-axis 
a 

direction 

v    Component  of  equilibrium  ship  speed  in  y-axis  direction 
e 

X    Hydrodynamic  and  aerodynamic  force  component  in  x-axis 
direction 

X    Aerodynamic  force  component  in  x-axis  direction 
a 

X    Hydrodynamic  force  component  in  x-axis  direction  due  to 

propeller 
X»    Derivative  of  hydrodynamic  force  component  in  x-axis 

direction,  with  respect  to  surge  acceleration 
X    Second  derivative  of  hydrodynamic  force  component  in 

x-axis  direction,  with  respect  to  sideslip  velocity  and 

yaw  angular  velocity 

X    Drag  coefficient 
o 

X    Second  derivative  of  hydrodynamic  force  component  in 
x-axis  direction,  with  respect  to  rudder  angle 

Y  Hydrodynamic  and  aerodynamic  force  component  in  y-axis 
direction 

Y  Aerodynamic  force  component  in  y-axis  direction 
a 

Y  Derivative  of  hydrodynamic  force  component  in  y-axis, 
with  respect  to  yaw  rate 

Y  Derivative  of  hydrodynamic  force  component  in  y-axis 
direction,  with  respect  to  sideslip  acceleration 

Y»    Derivative  of  hydrodynamic  force  component  in  y-axis 
direction,  with  respect  to  rudder  angle 


S  Rudder  angle 

S  Equilibrium  rudder  angle 

e  Density  of  water 

e  Density  of  air 
a 

ty  Heading  angle  of  ship 

ty  Direction  of  wind  velocity  relative  to  ship 

\p  Direction  of  wind  velocity  relative  to  earth 

ft  Angular  velocity  of  the  body  about  the  origin 
ft.=  ip  +  jq  +  kr 


10 


DIMENSIONLESS  FORMS 


Quantities 


Typical  Symbol 


Typical 
Dimensionless  Form 


Length 

Force 

Moment 

Mass 

Angular  velocity 

Static  force  rate 

Static  moment  rate 

Rudder  force  rate 

Damping  force  rate 

Damping  moment  rate 

Inertial  coefficient 

Moment  of  inertia 

Acceleration 

Velocity 

Angular  acceleration 

Density  of  air 

Inertial  coefficient 


Y 

N 
m 
r 
y 

N 


N 


u 
u 

r 

c 

N- 


X'    = 

r 

y  = 

N'    = 
m'    = 

v  •        — 


Y'    = 

v 

N'    = 
v 

y    = 

s 

Y'    = 
r 

»•  = 

i 
Y«    = 


I' 
z 

u' 
u' 
r' 

«; 

i 


xr/£ 

Y/^2U2 
N/^3U3 

r£/U 

V2^2U 
Nv/^3U 

Yg/^U2 

Y   /^3U 
r    2 

N   /^U 
r    2 

uJl/U2 

u/U 

r£2/U2 


c?a/« 
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I.   INTRODUCTION 

In  the  last  twenty  years  many  mathematical  derivations 
have  been  performed  in  the  area  of  motion  of  vehicles  in  a 
fluid  medium.  Some  of  these  mathematical  models  have  been 
derived  using  the  six  degree  of  freedom  equations  of  motion. 

In  this  study  of  six  degree  freedom  equations  are  re- 
duced to  three  degrees  of  freedom  with  several  assumptions. 
Low  ship  speed  is  assumed  and  coupling  effects  of  roll, 
pitch  and  heave  motion  on  to  the  horizontal  plane  motion 
are  assumed  negligible.   It  is  also  assumed  that  atmospheric 
motion  is  uniform. 

Aspects  of  the  problem  are  simultaneous  consideration  of 
hydrodynamic  and  aerodynamic  effects  as  well  as  nonlinearity 
in  the  external  disturbances  (especially  in  aerodynamic  terms) . 
Equilibrium  conditions  and  perturbation  equations  are  derived 
from  the  mathematical  model.   Nonlinear  parts  of  the  mathe- 
matical model  are  linearized  with  respect  to  equilibrium 
conditions . 

The  effects  of  wind  disturbances  on  ship  motions  are 
studied.   A  digital  computer  (IBM  model  360/67)  was  used  to 
solve  the  transient  response  in  yaw  and  sway  for  three  dif- 
ferent wind  directions  (head,  beam,  stern)  and  constant  wind 
velocities  with  constant  ship  speed.   Also  the  eigen  value 
theorem  was  used  to  study  ship  stability.   All  these  com- 
putations were  made  for  the  unsteered  (zero  degree  rudder 
deflection)  and  steered  case  (with  manual  as  well  as  auto- 
matic control) . 
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When  wind  velocities  are  several  times  the  ship  speed, 
even  moderate  winds  create  difficulty  in  controlling  the 
ship  course  and  stability  (at  low  ship  speeds) .   These  ef- 
fects may  or  may  not  be  controllable,  since  the  hull  hydro- 
dynamic  reactions  are  reduced  and  ship  motion  is  greatly 
effected  by  wind  force,  propeller  thrust  and  coupled  rudder 
force.   This  situation  can  cause  unusually  large  hull  drift 
angles  extending  to  between  plus  or  minus  180  degrees.   The 
usual  hull  hydrodynamic  force  representation  is  not  applica- 
ble beyond  15  to  20  degrees. 
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II.   THE  MATHEMATICAL  MODEL 

A.   COORDINATE  SYSTEMS  AND  SHIP  MOTION 

The  coordinate  system  is  defined  with  respect  to  the 
ship  itself.  It's  coordinate  origin  is  chosen  to  be  the 
center  of  gravity  of  the  ship  as  shown  in  Figure  1. 

Positive  X  axis  (lengitudinal  axis)  is  directed  to  the 
forward  direction  along  the  center  line  of  the  ship. 

Positive  Y  axis  (transverse  axis)  is  directed  to  star- 
board, perpendicular  to  the  plane  of  symmetry. 

Positive  Z  axis  is  directed  downward,  and  lies  in  the 
vertical  symmetry  plane,  perpendicular  to  the  X-Y  plane. 

Yaw  is  a  retatienal  motion  of  the  ship  about  the  Z  axis. 
(Assume  that  the  ship  is  moving  in  the  positive  X  direction 
and  that  a  moment  is  applied  about  the  Z  axis  tending  to 
rotate  or  turn  the  ship  in  clockwise  direction  as  indicated 
in  Figure  2  which  shows  the  positive  yaw) . 

1.   Equations  of  Motion  (on  the  horizontal  plane) 

The  dynamic  response  in  six  degrees  freedom  for  the 
motion  of  a  rigid  body  [7]  is  given  by  equation  (1) . 

X  =  W[U-RV+QW-X^  (R2+Q2)+Y/r  (PQ-R)+Z^  (PR+Q)  ] 

D  DO 

Y  =  W[V-PW+RU+Xr (R+PW)-Y. (P2+R2)+Z^ (RQ-P) ] 

b  fa  fa 

Z  =  W[W-QU+PV+X, (PR-Q)+Y. (P+QR)-Z^ (Q2+P2) ] 

b  b  b 

L  =  PIx+(Iz-Iy)QR+PI[Y6(W-QV+PV)-Z6(V-PW+RU)  ] 
M  =  QIY+(IX-IZ)PR+M[Z6(U-RV-QW)-X6(W-QU+PV) ] 
N  =  RIz+(Iy-Ix)PQ+pj[X6(V-PW+RU)-Y6(U-RV+QW)]  (1) 
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These  equations  represent  the  reaction  of  the  rigid 
body.   They  do  not  include  external  forces  such  as  rudder 
force,  wind  force,  aerodynamic  and  hydrodynamic  force  of 
water  etc.   Applied  external  hydrodynamic  and  aerodynamic 
forces  are  considered  to  be  functions  of  hull  motion,  rudder 
angle,  propeller  and  wind. 

Assuming  that  roll,  pitch,  and  heave  motions  may  be 
neglected  on  the  horizontal  plane  and  that  xc ,  y  c ,  zr    are 

boo 

zero  for.  the  center  of  gravity  at  the  origin  of  the  system 

.  ■*■"?■    -f 

and  noting  that  on  the  horizontal  plane  U  =  iu  +  jv    (linear 

velocity) ,  then  equation  (1)  becomes 

X  =  jn[{j  "  RVl  (Surge) 

Y  =  w[v  +  RUl  (Sway) 
Z  =  0 

L  =  0 
M  =  0 
N  =  RIZ  (Yaw)  (2) 

These  equations  are  nonlinear.   Analytic  solution  and 
computer  simulation  are  very  difficult.   But  all  nonlinear 
terms  might  be  linearized  by  using  the  Taylor  series  expansion 

Briefly,  small  perturbations  are  considered  for: 

X  =  m(u  -  rv) 

Y  =  m(v  +  ru) 

N  =  rl  (3) 
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The  external  force  and  moments  to  be  included  are:  [4] 

X  =  X(hull)  +  X (rudder  and  prop.)  +  X(wind) 

Y  =  Y(hull)  +  Y  (rudder  and  prop.)  +  Y(wind) 

N  =  N(hull)  +  N (rudder  and  prop.)  +  N(wind)       (4) 

So  N,  X,  Y  represent  total  hydrodynamic  and  aerodynamic  forces 
and  moments  which  are  generated  by  the  ship  motion,  rudder, 
propeller  and  wind.   The  following  assumptions  permit  elimina- 
tion of  relatively  unimportant  terms  in  the  hydrodynamic  and 
aerodynamic  expressions:  [1] 

1.  Ship  asymmetry  with  respect  to  the  centerline  plane 
is  minor  (this  permits  elimination  of  odd  ordered  terms  for 
the  X  force  expansion  and  even  ordered  terms  from  the  Y  and 

N  expansion) . 

2.  All  terms  higher  than  third  order  are  unimportant. 

3.  Rudder  force  and  moment  derivatives  of  higher  than 
first  order,  and  effects  of  rudder  angular  rate  are  negligible. 

4.  The  lengitudinal  component  of  ship  speed  is  much 
greater  than  the  lateral  component  (u>>v) . 

5.  Based  on  experimental  and  analytical  results,  non- 
linear hydrodynamic  terms  in  the  yaw  and  sway  equation  for  the 
stability  analysis  are  considered  minor,  except  for  certain 
nonlinear  hydrodynamic  terms  in  the  surge  equations  such  as 


X   ' 
vr 


6.   Atmospheric  motion  is  uniform. 

Thus  equation  (4)  may  be  written  as:  [1] 

2  2    2 

X  =    C,  u+C.u   +C0vr+C.u    5   +x,v  U   +x 
12  3  4  laap 
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Y  =  B,v+B„uv+B0ur+b  U26+Y.v  U 
^.2     3    y     1  a  a 

N  =  U+Luv+Lur+A  v  {+N.U  v  +N  v  U  (5) 

All  hydrodynamic  and  aerodynamic  constants  are  defined  in 
the  appendis. 

The  hydrodynamic  effects,  forces  are  defined  as  a 
function  of  the  acceleration  and  velocities  of  the  hull. 

The  rudder  deflection  will  also  produce  an  additional 

2 

resistance,  an  axial  force  proportional  to  S  .   This  force  is 

expected  to  be  small.   The  lateral  force  and  moment  are  pro- 
protional  to  S.  [5] 

Forces  and  moments  due  to  wind  acting  on  the  above 
water  protions  of  the  hull  act  on  the  ship  in  the  horizontal 
plane.   It  is  necessary  to  determine  the  relative  wind  speed 
and  direction  with  respect  to  the  ship  in  it's  own  reference 
frame . 

Magnitude  and  direction  of  wind  is  shown  on  Figure  1. 
The  relative  wind  speed  components  along  the  ship  axes  are 

u  =  u+U   cos  (iji  +i|»)       in  X  direction 
a      a       a 


v  =  v_u   sin  (d,  +*)       in  Y  direction 
a      a       a 


Then 


Ua  =  (ua  +  vl)k    '   *a  =  t™1—  <6> 

ua 

U   and  i>      are  wind  velocity  and  direction  along  the  earth 
a      a 

axis. 
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Sk.p    rtfe-rence.  frame. 


V 


On'oin    Of   tnfc 


>>i 


Figure   1.      Relation  of  Wind  Velocity   and  Direction   to   the 
Ship  Axes. 
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From  equations  (3) ,  (4) ,  and  (5) ,  accelerations  in 
yaw,  sway,  and  surge  are  determined: 
So 

X  =  u(m-xu)  -  mrv 


Y  =  v(m-yv)  +  mru 


N 


r(I   -  N-) 
z    r' 


(7) 


The  relationship  between  velocities  of  the  ship  with 
respect  to  the  earth  axis  and  velocities  of  the  ship  with 
respect  to  the  body  axis  is:  [3] 

Y  (t)  =  u  sin  (ij;)  +  v  cos  (^) 


also 


X  (t)  =  u  cos  (i|j)  -  v  sin  (ip) 


<Mt)  =  r 


(8) 


Where  ^  is  the  heading  angle  of  the  ship  relative  to  the 
original  path. 


oriQinoA        pQcth 


Actual     PCkxVi 


Figure  2.   Positive  Direction  of  Yaw 
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Y  =  /(u.sin  U)    +  v.cos  (40  )  dt  +  Y  (o) 

XQ  =  /(u.cos  {ty)    -   v  sin  (40)  dt  +  Xq(o)  (9) 

when  \p    is  small: 

sin  (ri>)-i>    ,        cos  (4>)~1    U^Ue  (equilibrium  case)    (10) 
So 

Y  (t)  =  Ue  •  i>   +   v 
o  Y 

X  (t)  =  Ue  -  v.i>  (11) 

Combining  equation  (5)  and  (7)  and  putting  in  nondimensional 
form.  [2] 

£il3  (m'-xu1  )YI=%12  [  (x'+x   '52)u2+Mx   ,+m,)vr+p'x,U2+x,U2] 
2  2     o   ss    '       vr        Ka  a  a   p 

^3  (m'-y  '  )v=^Jl2  [y'uv+Jl  (y'-m' ) ur+y ' 6u2+e 'y 'U2 ] 
z  v  2  v  r  aaa 

£jt5  (I'-N«  '  )r=;k3  [N'uv+JlN,ur+Nl6u2+elN,U2  (12) 

2  2     r   '      2  v  r  aaa  x      ' 


where 


r%3U2N'=N    =^3U2N' 
2  a   a      a    2  a 


N1  N1 

N:=p'U'2(0aisin(2^a)+9a2sin(^a)) 
a      a   a      z  a      z  a 


ail2U2y,=y  =^2U2y' 
z  aJa  Ja    2  -'a 


y'=p»U,2Y'    sin(^    ) 
Ja   Fa  a     ai  ra 


^2U2X'=X    =£il2U2X' 
2  a  a      a   2  a 


x'=P;u;2x:  cos(4)  (13) 

a      a   a      a  2  a 
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p»  =  p  /p   ,    u-U   since  u>>v) 

3.  a 

Expanding  in  Taylor's  series  around  the  equilibrium  conditions 

u  =  ue  v  =  ve  6  =  6e 

r  =  x  =  x  =  u  =  v=   =0  (14) 

Substitute  equation  (14)  into  equation  (12) 

N'ueve  +  N'6e  ue2  +  N'ue2  =  0 
v        6  a 

y'ueve  +  yl6e  ue2  +  y'  ue2  =  0 

Jv  •      2 6         J a 

x'ue2  +  x   '6e2  ue2  +  x*  ue2  +  x'ue2  =  0  (15) 

o       66  a        p 

Sideslip  velocity  and  rudder  angle  in  the  equilibrium  condi- 
tion can  be  expressed  as:  [1] 


Y  'N  '  -  N  'Y  ' 
..  ,     a   6     a   6 
Ve'  =  


N  'Y  '  -  Y  'N  ' 

v   6     v   5 


N  'Y  '  -  Y  'N  ' 
6e  =  -2—Z <L_Y_  (16) 

N  'Y  *  -  Y  'N  ' 

v   6      v   6 

Perturbation  terms  (small  deviation  from  the  initial  straight 
course)  may  be  denoted  by  bars.   So  actual  motion  is 

ue  =  ue  +  U     V  =  ve  +  V     S  =  Se  +  S     ^  =  ? 

i  -  i  u  =  u       v  =  v  (17) 


Perturbations  in  aerodynamic  quantities  due  to  $,  v'  and  u* 
are: 
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3N'               »N' 

N'    =         a    J    +   —2.  v' 

+ 

3N1 
a 

u' 

a         d4>                3V' 

3u' 

3Y'              3Y' 

Y'    -         a    J    +         a   v' 

+ 

3Y' 
a 

u1 

a         8  4,                3v' 

3u' 

ax'         ax' 

3X' 

X'    =   — -  5   +   — -  v'    +   — -  U'  (18) 

3ip  3v'  3u' 

Where 

3N' 

— -  =   p '  [  (N '    cos    2^      +  N *      cos   $    )     (U ' 2    -   u ' ) 

3.        3.  i  a  a  2  a  a  a 


34> 


-2U1    sin   ty       (— x    sin   2ip      +  N'      sin   ti)  ] 
a  a      ^  a  a  o  a 


3N'  N' 

— -  -   -pi   U'[(N'      cos    2*      +  N'    cos   iL    )    +   2  (— !    sin   24, 
~     i  a      a        a^  a  a2  a  ~  a 


+  N'      sin   \b    )    sin   \b    ] 
a2  3  a 


3N' 

— -  =   pl    U1  [  (N1       cos    2i>      +   N'       cos    ty    )  (-sin   i>    ) 
i    i  a      a        a^  a  a2  a  a 

N' 

+  2  (— l    sin  2^   +  N'   sin  i>    )    cos  if*  ] 
~         a    a2      a       a 


3Y' 

-  Pl  yl  [cos  ij>  (U'2-ul)  -  2U'  sin2  *  ] 

j\       a  a  i  a   a    a      a       a 


3Y' 

—2-  =  -p'  Ul  Yl   (1  +  sin2  *  ) 
9v'      a   a   a2  a 


3Y.' 

gui    "~  a  ~ a  ~ a i     T  a     Ta 


=  p*  U'  Y"   cos  ty      sin  i/>, 
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3X' 

— -  =  -  p'  X'   sin  ty   (U'2  +  u') 
.       Ka   a,      ra   a     a7 

3x  l 


— -  =  -p '  X'   u1  cos  i>      sin  4> 
3v«      a   aj   a      a      a 


DX ' 

-  =  p'  X'   U'   (1  +  COS2  $    )  (19) 

8u'     a   aj   a  a 

The  coefficients  in  equation  (19)  where  taken  from  reference 

[2],   Aerodynamic  coefficients  are  shown  in  Figure  3. 

Substitute  equation  (17)  and  (18)  into  equation  (12) 

and  obtain  the  perturbation  equation  in  the  following  form: 

2x'   2x'  2X ' 

(m'-x-  '  )u'  =  (5   ,+m,)v  'r'  + — -5+ — -v'  +  (2X'+2XJC  !<S2  + — -) 

2x    2V       °    66  e  2u' 


u,+2X.,'6  6 
66   e 

2Y'       2Y1  2Y' 

(m'-y  'Jv'^y  l-m')r'+ — =+  (y  '  + — -)  v'+  (y  '  v'+2y  '6e+ — -) 
V        2x    v   2V       V  e    6    2u' 


u '  +v  J  6 


*6 


2N'       2N'  2N' 

(I'-N-  ' )r'=N'r'+ — -J+(N'+ — S.) v'  +  (N1  v'  +2NJ6  + — *) 
Z   r      r    2*      V  2v'       v  e    6  e  2u' 


u'+Njfi  *  (20) 

The  mathematical  model  which  is  shown  in  the  above  equations 
is  developed  only  for  the  forward  ship  motion.   The  stern 
ship  motion  equations  are  written  by  putting  negative  signs 
in  front  of  some  of  the  coefficients. 
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Example:   N   (derivative  of  hydrodynamic  yaw  moment  with 
respect  to  sideslip  velocity)  takes  a  negative  value  for  the 
forward  motion,  positive  value  for  the  stern  motion. 

Y   (derivative  of  hydrodynamic  force  component  in 
Y-axis,  with  respect  to  yaw  rate)  takes  negative  value  for 
the  forward  motion,  positive  for  the  stern  motion  but  is  not 
a  sensitive  parameter  since  m'  is  the  predominant  term  in 
(Yr'  -  m'). 
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III.   CONTROL  STUDIES 

The  linearized  equation  of  motion  was  derived  and  found 
in  equation  (20) .   In  this  part  two  different  problems  were 
simulated.   One  of  them  was  manually  steered  ship  motion  in 
wind  disturbance.   The  other  was  automatically  steered  ship 
motion  under  the  same  effect.   So  computer  program  (II)  and 
(III)  were  developed  for  these  purposes  respectively.   The 
computer  model  is  shown  in  Figure  4,  for  these  programs. 
Program  (II)  and  (III)  are  similar  to  each  other.   There  is 
only  one  difference  in  section  4  of  both  programs.   Program 
(II)  was  integrated  using  the  second  order  equation,  however 
program  (III)  was  integrated  with  the  third  order  equation. 
This  difference  was  caused  by  the  feedback  loop  (page  38  ) . 
Also  the  Laplace  Transform  Technique  could  have  been  used  to 
obtain  the  system  description  for  simplicity  and  system 
characteristic  equation  for  the  parameter  plane  program.   The 
characteristic  equation  was  calculated  in  the  computer  pro- 
gram.  Also  three  more  additional  and  different  mathematical 
sections  were  calcualted  in  the  same  program.   One  of  these 
sections  included  the  parameter  calculation  from  equation 
(6)  through  equation  (20) .   Secondly  the  transfer  functions 
of  yaw  and  sway  were  calculated  from  the  three  degree  of 
freedom  equation.   The  last  section  was  the  polynomial  cal- 
culation for  the  eigen  values. 

The  last  program,  (parameter  plane  (IV) ,  was  used  to  get 
optimal  values  for  yaw  gain  constant  and  yaw  rate  gain  con- 
stant using  the  constant  zeta  (function  omega)  and  constant 
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omega  (function  zeta)  curves  for  automatically  steered  ship 
motion.   According  to  reference  [1]  and  [6] ,  these  optimal 
values  were  chosen  from  Figures  40  through  44  assuming 
zeta  =  6,  omega  =  11.5  for  varying  wind  direction  and 
velocities. 

The  transform  of  the  linearized  equations  of  motion  (20) 
from  the  T  domain  to  the  S  domain,  for  all  initial  conditions 
equal  to  zero  are: 

(ABC.£2-AAD.(S-BFA)  i>-  (BGD.tf)—  -(BGC.£)~  =  AAE  .  6 

6  6 


(-BGE.S-BFD)J- (BGG.S)-  + (ABD.£2-BGF .£) -  =  AAC. 6 

6  6 


(21) 


(-CAA.*-BFG)i£+(ABE..S2-CAB.S)~  -(BGA.3)-  =  CAC. 6 


In  equation  (21)  different  letters  are  used  for  simplicity. 

For  example:    ABC  =  I'  -  N- *  ,    AAD  =  N'   ,    BFA  =  8Na  , 

z     r  r  3i|> 


ABE  =  m1  -  x*  *   ,   etc. 

u    ' 


let 


H  =  c 
s   C   ' 


1=  B 
S 


Applying  Cramer's  rule,  the  transfer  function  for  —  and  — 


are: 

s 


AAE 
AAC 
CAC 


-BGD . S  AAE 

-BGG . S  AAC 

ABE.S2-CAB.S    CAC 


(ABC.S2-AAD.S-BFA)    -BGD.S  -BGC . S 

(-BGE.S-BFD)  -BGG.S  ABD.S2-BGF.£ 

(-CAA.S-BFG)  ABE.S2-CAB.S    -BGA.s      I 


(22) 
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The  open  loop  transfer  function  is: 


$         WA.S2+WC.S+WF 

6    SSA.Sit+SSB.S3+SSE.S2+SBS.S+SES 


(XI) 


Where 


WA  =  -  ABC.ABD.ABE    etc. 

All  terms  were  calculated  using  program  1. 

In  this  case,  the  results  are  a  second  order  numerator 

and  fourth  order  denominator. 

B(s) 


Transfer  function 


S(s) 


,  was  also  obtained  using  Cramer's 


rule. 


B 


ABC.S2-AAD.S-BFA 

-BGD.S 

AAE 

-BGE.S-BFD 

-BGG.S 

AAC 

-CAA.S-BFG 

ABE. S2 -CAB. S 

CAC 

ABC.S2-AAD.S-BFA 

-BGD.S 

-BGC.S 

-BGE.S-BFD 

-BGG.S 

ABD.S2-BGF.S 

-CAA.S-BFG 

ABE.S2-CAB.S 

-BGA.S 

(23) 


The  transfer  function  between  B  and  £  is: 


B        WWA.S3+WWC.S2+WWF.S+WVB 

6    SSA.S5+SSB.St++SSE.S3+SBS.S2+SES.S 


(24) 


where 

WWA  =  -ABC. AAC. ABE   etc. 
Again  all  terms  were  calculated  using  program  1 
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For  the  lateral  displacement,  sway  is  defined  by 


S  _ 


WWA. S  3+WWC . S2+WWF . S+WVB 
SSA.S5+SSB.St++SSE.S3  +  SBS.S2+SES.S 


WA .  S  2  +WC  .  S+WF 


SSA.S4+SSB.S3+SSE.S2+SBS.S+SES 


s 


f 


y 


s 


v 


The   transfer   function   is 

Y        WWA.S    +(WWC+WA)S    +  (WWF+WC)  S+WVB+WF 

S         SSA.S    +SSB.S    +(SSE+SSA)S    + (SBS+SSB) S    +(SES+SSE)S    +SBS.S+SES.S 

(25) 

A.   MANUALLY  STEERED  SHIP  MOTION 

As  mentioned  before,  the  goal  of  this  section  is  to  get 
the  ship  trajectories  in  constant  wind  effect  with  low  ship 
speed.   For  this  purpose  computer  program  (II)  was  developed 
and  applied  to  equation  (20).   Two  situations  were  considered: 

Situation  1:   zero  degree  rudder  deflection  (unsteered) . 

Situation  2:   positive  15  degree  rudder  deflection. 
All  these  situation  were  studied  assuming  that  the  ship  moved 
in  the  forward  direction  at  a  1  mile/hr  speed  in . a  calm  sea. 
Each  of  the  wind  velocities  was  applied  suddenly. 

Computer  program  II  responses  are  shown  in  Figures  5 
through  30  for  yaw  and  sway.   These  results  were  also  tabu- 
lated in  Tables  1  through  4.   The  first  test  on  the  unsteered 
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situation  is  shown  in  the  computer  outputs  from  Figure  5  to 
Figure  24  for  three  different  wind  directions  and  various 
velocities.   Tabulated  values  are  shown  on  Table  1  and  2, 
for  yaw  and  sway. 

The  effects  are  examined  by  referring  to  the  tables. 

Casel:   Head  Wind 

As  is  shown  in  Tables  1  and  2  and  Figures  5-10,  effects 
of  the  head  wind  were  small  on  the  straight  forward  ship 

motion.   In  this  test,  a  5  mile/hr  wind  velocity  is  applied 

-3 

to  the  ship  which  created  -10    rad.  variation  on  the  course 

_3 
and  drift  was  -7.4x10    ft.  in  5  minuted. 

Case  2 :   Beam  Wind 

Secondly,  beam  wind  disturbances  were  examined  for  various 
velocities.   According  to  Tables  1  and  2  and  Figures  11  -16, 
the  ship  greatly  deflected  from  a  straight  course  where  the 
wind  velocity  exceeded  t  miles/hr.   Deviation  was  positive 
22  degree  (.395  rad)  on  the  straight  course  and  positive 
9.2  ft  drift  in  5  minutes  when  the  7  miles  wind  velocity  was 
applied,  as  shown  in  Figures  15  and  16. 

Case  3:   Stern  Wind 

In  this  part,  the  ship  was  disturbed  with  stern  wind 
which  created  great  difficulty  on  the  straight  course.   Five 
miles/hr  wind  velocity  created  -2.29  degree  (-.04  rad)  devia- 
tion and  r.28  feet  drift  in  5  minute.   But  7  miles/hr  wind 
velocity  causes  the  ship  to  deviate  by  -332.3  degree  from 
the  straight  course  in  4  minute  10  second.   Drift  was  -110  ft 
in  4  minute  but  this  value  was  reduced  to  -15  ft  in  5  minute. 
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Briefly  the  system  was  unstable  and  the  ship  was  turned 
around.   There  is  no  chance  to  control  the  ship.   Drift  was 
negative  110  ft  in  4  minute  but  this  value  was  reduced  the 
negative  15  ft  in  5  minute.   The  summarization  of  this  test 
can  be  stated  as  the  following:   computer  program  (II)  re- 
sponses were  tabulated  for  3  different  wind  directions  and 
various  velocities.   According  to  the  table,  wind  distur- 
bance was  not  so  important  on  the  course  from  any  direction 
until  5  miles/hr  wind  velocity  was  reached.   Critical  points 
were  5  and  7  miles/hr  wind  velocities  for  stern  and  beam 
wind  respectively. 

Second  test:   15  degree  rudder  displacement  was  applied 
under  the  beam  and  stern  wind  disturbances.   The  main  pur- 
pose was  to  show  the  effect  of  the  rudder  displacement  on  the 
ship  course  with  wind.   As  is  known,  positive  rudder  displace- 
ment causes  negative  yaw.   All  responses  are  shown  in  Figures 
25-30  and  tabulated  values  are  shown  in  Tables  3  and  4.   Ac- 
cording to  the  Tables  3  and  4 ,  when  the  ship  is  effected  5 
miles  beam  wind  velocity  and  positive  15  degree  rudder  dis- 
placement.  So,  ship  was  deflected  -327  degree  from  the 
original  path  in  5  minute.   Drift  was  80  ft  in  2.5  minute 

and  -100  ft  in  5  minute.   As  mentioned  earlier,  5  miles/hr 

-3 
beam  wind  with  zero  rudder  displacement  has  created  8  x  10 

rad.  deviation  from  the  original  path,  as  shown  in  Table  1. 

But  the  same  condition  with  15°  rudder  displacement  has  -327° 

deviation  from  the  original  path.   In  this  case,  the  ship 

was  so  effected  with  rudder  force  and  started  to  turn  around 

the  path. 
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Figure   WD. PI.  (degree)   WD.VEL.  (miles)   Time (sec)   Yaw(rad) 


180 

3 

300 

-4 

180 

5 

300 

-18 

90 

5 

300 

-5.7 

Table  3.   Variation  of  Yaw  for  Manually  Steered  Ship  (S=15°) 


Figure   WD. PI.  (degree)  WP.VEL.  (miles)   Time  (sec)   Sway (ft) 

180  3  <240      {?18 

180  5       nonlinear  oscillation 

qn  S  I150       f80 

90   .  5  {300       {-100 


Table  4.   Variation  of  Sway  for  Manually  Steered  Ship  (S=15°) 
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The  other  situation  is  a  stern  wind  disturbance  with  15° 
rudder  displacement.   Also  known  from  Table  1,  5  miles/hr 
stern  wind  velocity  with  zero  rudder  displacement  has  created 
-.04  rad.  deviation  from  the  original  path.   In  this  case, 
15  degree  rudder  displacement  created  -18  rad.  deviation  in 
5  minute.   Ship  was  turned  360°  around  the  course  in  2  minute 
and  drifted  with  divergent  oscillation. 

All  these  tests  were  made  with  ship  speed  being  1  mile/hr 
in  the  calm  water  and  values  of  wind  velocities  were  suddenly 
applied. 

B.   AUTOMATICALLY  STEERED  SHIP 

Earlier  tests  show  that  the  ship  was  effected  by  wind 
and  rudder. 

In  this  section,  another  test  was  developed  to  minimize 
the  deviation  from  the  original  path.   For  this  reason 
rudder  deflection  is  used  and  is  written  as  a  function  of  ty. 

The  equation  for  the  rudder  deflection  is: 


.6   =  a   5d  +   5d  =  k      ty   +  k\\> 


where 

a  =  rudder  response  time  constant  (time  lag) 

k  =  yaw  gain  constant 

k.  =  yaw  rate  gain  constant 


Thus,  the  system  would  be: 
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•<S4I 


PL* 


NT 


Pinnr 


KtS 


K 


These  control  effectors  are  the  function  of  some  parameter 
that  could  be  measured  by  sensors,  usually  a  motion  parameter 


*J*/^Y& 


i 


CX.S+' 


Let 


1 

—  =  a, 

a    1 


vva.  sl+  vvc  S  +■  wr 


SSO    SV+SS6  SV  SSf.sV  3/sS.^-f  ScS 


kt5    +  K 


£        .  WA.S2+   WC.S+WF 


a. 


5d         SSA.S"+SSB.S3+SSE.S2+SBS.S+SES 


S+a. 


H    = 


K    +    KtS 


(26) 
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So  the  characteristic  equation  is 

SSA.S5+(SSB+cx,/SSA)  .  S*  +  (SSE+a  SSB-a  .  WA.K  )  S  3 
+  (SBS+a;..SSE-a  .WC.K  -a  .WA.K)S2+(SES+a1.SBS 
-o,.WF.K   -a1.WC.K)S+(a    .SES-a    .WF.K)    =    0  (27) 

For  simplicity,  let  B  =SSA,   B  =SSB+a, «SSA  etc. 

B,  -S^B-S^B  -B..K  )S3+(B  -B  -K  -B_K)S2 
1      2      34   t       56t7 

+(B8-B9.Kt-B10.K)S+(B11-B12.K)  =  0  (28) 

This  characteristic  equation  was  calculated  using  computer 
program  1.   Computer  program  III  was  used  to  get  transient 
responses  in  yaw  and  sway  for  the  automatically  steered  ship. 
Each  curve  of  yaw  or  sway  was  developed  with  different  param- 
eter values  (for  K  and  K, ) ,  for  different  wind  speed  and 
direction. 

The  problem  is  to  minimize  the  variation  on  the  ship 
course.   For  this  reason,  the  parameter  plane  program  is  ap- 
plied to  the  system  to  determine  the  best  value  of  yaw  gain 
constant  and  yaw  rate  gain  constant. 

The  responses  of  these  are  shown  from  Figure  40  to 
Figure  44  for  several  situations.   The  values  chosen  are 
the  minimum  values  shown  on  the  curves  according  to  refer- 
ence [6]  . 

The  parameter  plane  outputs  are  developed  by  applying 
equation  (27)  to  program  IV.   Table  5  shows  the  best  values 
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of  the  yaw  gain  constant  and  yaw  rate  gain  constant  for  the 
automatically  steered  ship. 

All  responses  of  the  automatically  steered  ship  are 
shown  in  Figure  31  through  33  and  tabulated  on  Tables  6  and 
7. 

According  to  the  Tables  (6  and  7) ,  the  variations  on 
the  course  is  minimized  for  the  low  wind  velocities.   In 
the  head  wind  disturbance,  course  variations  were  reduced 
to  0  rad.  in  5  minute. 

One  mile/hr  wind  velocity  response  reached  0  rad.  in 

20  sec,  3  mile/hr  gave  0  rad.  response  in  65  sec.   But  5 

-4 
miles/hr  wind  velocity  response  was  1x10   rad.  in  5  minute 

For  the  beam  wind  effection,  3  miles/hr  wind  velocity 
was  applied  to  the  ship.   Deviation  was  0  rad.  in  70  sec. 

A  stern  wind  was  also  tried,  1  mile/hr  wind  velocity 
caused  a  deviation  which  was  corrected  to  0  rad.  in  65  sec. 

As  is  seen  above,  the  ship  trajectory  was  kept  on  the 
original  path.   This  desired  course  or  original  path  was 
assumed  to  be  due  North.   Also  the  same  test  condition  was 
applied. 
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Figure 

WD.DI. (de 
0 

;gree) 

WD 

.VEL. < 
1 

mi 

les) 

Time (sec) 
r  10 

Yav/  (rad) 
r-lxl0~4 

^300 

-0 

180 

1 

,  50 

^00 

,-lxlO"4 
l~0 

0 

3 

r   15 

r-lxlO"4 

l300 

l~o 

90 

3 

r   20 

^00 

r.85xl0"4 

l~o 

0 

5 

,   30 

,lxlo""4 

*300 

i. 5x10-4 

Table  6.   Variation  of  Yaw  for  Automatically  Steered  Ship. 


Figure   WD.DI.  (degree)   WD. VEL.  (miles)   Time (sec)   Sway (ft) 


0 

1 

300 

10xl0"4 

180 

1 

300 

13xl0"3 

0 

3 

300 

11.5x10" 

90 

3 

300 

-5xl0~3 

Table  7.   Variation  of  Sway  for  Automatically  Steered  Ship. 
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IV.   STABILITY  ANALYSIS 

Stability  may  be  defined  in  the  following  way.   Estab- 
lish an  equilibrium  situation,  disturb  the  equilibrium  con- 
dition with  arbitrary  disturbances.   If  the  ship  returns  or 
tends  to  return  to  the  original  straight-line  condition  of 
equilibrium  after  a  disturbances  of  the  smallest  amount 
(infinitesimal  disturbance),  it  is  stable.   If  it  departs 
or  has  the  tendency  to  depart  from  the  original  straight- 
line  equilibrium  condition,  it  is  unstable. 

The  unstable  ship  will  go  in  to  a  starboard  or  port 
turn. 
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Stability  analysis  is  performed  by  solving  for  the  eigen 
values  of  the  system  matrix  (4x4)  which  is  defined  on  page 
46.   Reference  [3]  contains  straight  forward  algebraic  solu- 
tion of  the  linearized  equations.   According  to  reference 
[3]  a  complete  algebraic  solution  of  linearized  equation  (20) 
is  obtained  with  a  differential  operator.   Solution  result 
of  any  independent  variable  X.  is  defined  as: 

0,t       oyt  o   t 

X.,.v  =  C..e  L      +  C0.e  £      + +  C  .e  m 

J  (t)     13        2j  mj 

m        c,  t 

xj(t)  =  ^=1cije  (29> 

Where  C. .  are  arbitrary  constants  depending  upon  the  initial 

conditions  and  a,,  ff_,  ....  a   are  the  roots  of  the  deter- 

12        m 

minant  of  the  coefficients  in  the  linearized  equations  of 

motion. 

In  general  a,,  a0,  ....  a   are  different  in  value,  the 
3        1    2        m 

°lt  CJ2t  Omt 

terms  C,e    .  C~e    ,    ....  C  e  m  can  not  negate  one  another. 

12  m  J 

The  only  way  that  each  term  can  go  to  zero  with  increasing 
time  is  for  each  of  the  exponents  to  be  negative.   If  a  is 
a  complex  number  in  the  form  a=a+ib,  the  following  relation- 
ship hold  e   =e       =e   (cos  bt+i  sin  bt)  and  the  condition 
for  stability  requires  that  the  real  parts  of  a,,  cr_i  ....  a 
be  negative  as  they  are  complex  numbers  (the  imaginary  a, 
parts  of  number  indicate  the  angular  frequency  of  oscillation) 
The  Routh  oriterian  can  also  be  applied  to  determine  system 
stability.   If  we  look  at  equation  (20)  again,  the  automat- 
ically steered  ship  equation  would  be: 
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ABC.^-AAD.^-BFA.^-BGC.V-BGD.U  =  AAE.S 
-BGE.if-BFDij;  +ABD.V-BGF  V-BGG  U  =  AAC .  S 
-CAA.ijJ-BFGiji-  BGA.V+ABE.U-CAB.U  =  CAC.S 


QOQ.ip+QQO.^  =  QQQ.S  +  S 
Where 


QQO  =  kt, 


QOQ  =  k, 


QQQ  =  tr  =  a 


Equation  (30)  will  be  written  in  the  following  way 


ABC 

-AAD 

-BFA 

0 

-BGC 

0 

-BGD 

0 

-AAE 

0 

$+ 

-BGE 

• 

-BFD 

fr 

ABD 

• 

V+ 

-BGF 

V+ 

0 

• 

U+ 

-BGG 

U+ 

0 

• 

-AAC 

0 

-CAA 

-BFG 

0 

-BGA 

ABE 

-CAB 

0 

-CAC 

0 

QQO 

QOQ 

0 

0 

0 

0 

-QQQ 

-1 

(30) 


g=o 


(31) 

As  is  seen  from  equation  (31)  there  are  four  independent 

variables  namely  yaw,  sway,  surge  and  the  rudder  angle.   The 

form  used  in  equation  (29)  is  applied. 

m       ,     m       .     m       .     m       . 
ip=Z   if.e01  ,  V=E   V.eair,  U=Z   U.eai  ,  S=Z   A.e0!11   (32) 


i-1  * 


i=l 


i=l 


i=l 


Where  ty . ,    V.  ,  U.  and  A.  are  arbitrary  constants. 

Substitution  of  equation  (32)  into  equation  (31) : 


ABCC    -AAD. Q. -BFA 
1  1 


-BGEO . -BFO 

1 


-CAAQ . -BFG 


QQOai+QOQ 


V 


-BGC 


ABDO . -BGF 

1 


-BGA 
0 


V.+ 

1 


-BGD 
-BGG 


ABEQ . -CAB 

1 


U.+ 

1 


-AAE 
-AAC 
-CAC 


-OOQc^-l 


A.=  0 

x 


(33) 
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Solutions  of  equation  (33)  are  possible  for  only  values  of 
a.  which  satisfy  equation  (34). 

(ABC0.2-AADa.-BFA)  (-BGC)        (-BGD)        (-AAE) 


(-BGEa.-BFD)         (ABDa.-BGF)  (-BGG)        (-AAC) 
(-CAAa.-BFG)         (-BGA)        (ABEa.-CAB)  (-CAC) 


=  0 


(QQOO..+QOQ)  0  0        (-QQQc^-1) 

(34) 

Program  I  gives  the  solution  of  this  equation  as  fifth 

i 

order  polynomial.   For  different  situations  such  as  wind 

velocity  and  direction  other  polynomials  were  also  found. 

As  is  pointed  out  in  page    ,  it  is  important  for  the  real 

root  to  be  negative  for  system  stability. 

This  test  was  examined  in  two  situations: 

Situation  1:  constant  wind  directions  with  varied 
velocities . 

Situation  2:  constant  wind  velocities  with  varied 
direction. 

The  first  situation  was  divided  into  3  sections: 

Section  a:  constant  head  wind  with  several  velocities. 
Section  b:  constant  beam  wind  with  several  velocities. 
Section  c:  constant  stern  wind  with  several  velocities. 

The  second  situation  was  divided  into  2  section: 

Section  a:  constant  7-  miles/hr  wind  velocity  with 
several  directions. 

Section  b:  constant  9  miles/hr  wind  velocity  with 
several  directions. 

In  order  to  examine  the  stability  of  an  automatically  steered 

ship  in  wind,  all  these  situations  are  solved  with  computer 
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program  I,  tabulated  on  Tables  8-12  and  drawn  on  the 
CALCOMP  plotters,  from  Figures  4  5  through  49. 

Situation  1:   This  test  was  developed,  constant  wind 
direction  disturbed  the  ship.   Wind  velocities  changed 
each  time. 

a)  Head  wind  disturbance: 

As  is  shown  in  Table  8  and  Figure  45,  the  ship 
would  be  stable  in  varied  head  wind  velocities  until  13 
miles/hr  wind  velocity.   After  this  range,  the  ship  could 
be  unstable  on  the  horizontal  plane. 

b)  Beam  wind  disturbance: 

In  this  case,  the  stability  range  was  shorter 
(from  1  to  8  miles/hr) .   After  the  8  miles/hr  wind  velocities 
the  ship  could  be  strictly  unstable,  shown  on  Table  9  and 
Figure  46. 

c)  Stern  wind  disturbance:         , 

Stern  wind  effect  on  the  stability  might  be  the 
shortest  range,  as  shown  on  Figure  47  and  Table  10.   This 
range  included  from  1  to  6  miles/hr.   After  6  miles/hr  wind 
velocities,  the  ship  could  be  unstable. 

Situation  2:   In  this  case,  constant  wind  velocity  dis- 
turbed the  ship.   Wind  direction  was  changed  for  each  test. 

a)  In  7  miles/hr  wind  velocity  with  varied  direction. 
From  0  degree  through  180  degree  wind  directions  were  tried 
in  7  miles/hr  wind  velocity.   According  to  Table  11  and 
Figure  48,  the  ship  could  be  unstable  when  the  wind  dis- 
turbed the  ship  from  30  degree,  105  degree  and  180  degree 
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wind  direction.   In  the  other  wind  direction  and  7  miles/hr 
wind  velocity  the  ship  trajectory  could  be  stable. 

b)  In  9  miles/hr  wind  velocity  with  varied  wind 
direction.   As  is  shown  on  Figure  49  and  Table  12,  9  miles/hr 
wind  velocity  with  each  direction  were  applied  to  the  system. 
The  ship  could  be  stable  when  the  wind  direction  was  head, 
60  degree  -  75  degree  and  120  degree  -  135  degree.   For  the 
rest  of  the  directions  (0°  -  180°)  the  ship  trajectory  was 
unstable. 

All  of  these  results  show  that  some  of  the  wind  direc- 
tions and  velocities  played  a  very  active  role.   These  direc- 
tions are  30  degree,  105  degree  and  180  degree.   After  the 
9  miles/hr  wind  velocities,  many  of  the  wind  directions  had 
adverse  affects  on  stability. 
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WD.DI. 

Re  a 

Re  a2 

Re  a3 

Re  a. 

4 

Re  o5 

0° 

-7.21 

-.124 

-.014 

-11.11 

-60.25 

30° 

+3.46 

-.0104 

-.138 

-10.215 

-60.44 

45° 

-3.34 

-.082 

-.0114 

-16.42 

-60.38 

60° 

-4.9 

-.088 

-.035 

-10.63 

-60.3 

75° 

-4.92 

-.145 

-.0143 

-10.6 

-60.26 

90° 

■   -2.56 

-.148 

-7.8xl0~3 

-10.4 

-60.4 

105° 

+  3.16 

-.0419 

-.419 

-10.2 

-60.56 

120° 

-6.22 

-.0896 

-.0153 

-10.8 

-60.29 

135° 

-6.56 

-.118 

-.0249 

-10.92 

-60.28 

150° 

-.0113 

-.947 

-.198 

-10.331 

-60.364 

165° 

-1.6xl0"3 

-.956 

-.155 

-10.3 

-60.44 

180° 

+  .443 

-.0286 

-.0286 

-10.286 

-60.62 

Table  11.   Real  Roots  of  Characteristic 

Equation,  Wind  Velocity  =  7  miles. 
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WD.DI. 

Re  o 

Re  a2 

Re  a3 

Re  a. 
4 

Re  a5 

0° 

-5.08 

-.138 

-.024 

-10.68 

-60.28 

30° 

+8.53 

-.19 

-.0238 

-10.16 

-60.1 

45° 

+.1054 

-1.25 

-.0275 

-7.85 

-60.8 

60° 

-2.25 

-.073 

-.0732 

-10.4 

-60.36 

75° 

-2.14 

-.193 

-.0224 

-10.42 

-60.33 

90° 

+1.94 

-.027 

-.145 

-10.25 

-60.42 

105° 

+8.14 

-.056 

-.056 

-10.16 

-60.62 

120° 

-3.07 

-.07 

-.045 

-10.45 

-60.42 

135° 

-4.31 

-.132 

-.038 

-10.55 

-60.34 

150° 

+2.8 

-.0368 

-.17 

-10.24 

-60.4 

165° 

+4.1 

-.0128 

-.146 

-10.2 

-60.47 

180° 

+4.67 

-.054 

-.054 

-10.21 

-60.55 

Table  12.   Real  Roots  of  Characteristic 

Equation,  Wind  Velocity  =  9  miles. 
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V.   DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

The  procedures  for  system  course  control  and  stability 
analysis  on  the  horizontal  plane  were  described.   Applying 
the  datas  which  are  taken  from  some  references,  to  the 
three  degree  freedom  equation  gave  us  course  control  and 
stability  analysis  computer  outputs.   As  a  result  of  these 
outputs  some  comments  can  be  made.   First  of  all  the  un- 
steered  case  was  examined.   It  is  seen  that  the  stern  wind 
effects  the  straight  course  of  the  ship  the  most.   When  the 
wind  velocity  reaches  7  miles/hr  this  effect  becomes  ex- 
tremely high.   In  5  minutes,  the  deviation  of  the  course 
is  about  327  degree.   At  the  same  time  the  effect  of  the 
head  wind  on  the  straight  course  compared  to  the  stern  and 
beam  wind  seems  to  be  much  less;  which  is  seen  from  the 
computer  outputs . 

To  see  the  effect  of  the  rudder  deflection  on  the 
straight  course  there  should  be  positive  15  degree  applied 
to  the  rudder  for  beam  and  stern  wind.   Under  these  con- 
ditions negative. yaw  and  sway  were  found.   When  the  ship 
is  under  stern  wind,  applying  positive  15  degree  rudder 
deflection  will  create  extremely  negative  yaw  and  sway. 
As  it  is  seen  in  Figures  29  and  30,  when  the  stern  wind 
velocity  reaches  5  miles  in  50  seconds,  15  degree  rudder 
deflection  will  create  minus  145  degree  deviation  from 
the  course. 
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Computer  output  results  were  also  obtained  for  an  auto- 
matically steered  ship.   Parameter  plane  computer  programs 
were  applied  to  the  closed  loop  system  for  choosing  the 
best  value  of  yaw  rate  gain  constant  and  yaw  gain  constant. 
According  to  the  automatically  steered  ship  computer  out- 
puts, in  low  wind  speeds  the  course  variation  is  minimized. 
Sway  can't  be  controlled  very  well  but  the  value  of  the 
drift  has  decreased  compared  to  the  unsteered  case. 

In -a  final  study,  using  the  eigen  value  theorem,  system 
stability  and  the  critical  wind  velocity  are  found  and 
tabulated.   Computed  results  show  that  30  degree,  105  degree 
and  180  degree  wind  direction  are  important  on  the  ship 
course  in  X  miles/hr  wind  velocities.   Also  9  miles/hr  wind 
velocity  in  many  direction  were  so  effected. 

All  the  studies  done  above  are  made  for  a  ship  which 
has  a  speed  of  1  mile/hr  and  limited  special  full  loaded 
ship  form. 
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APPENDIX 

Hydrodynamic  and  aerodynamic  coefficients  where  used 
in  equation  (5) 

Al  =  I  *5  Nr'     Bl  =  !  **  V     Cl  =  !  *3  V 

A2  =  !  £3  V     B2  =  !  ll    V     C2  "  !  *2  Xo' 
A3  =  |  £*  Nr«     B3  -  |  £3  V     C3  =  |  ^  x^. 

A,  =  §  £3  N  '     B,  =  ^  I1    Y  ■     C„  =  %   £2  X   ' 
42      s       42      s       42      ss 

Where 

A.  =  Coefficients  of  hydrodynamic  moment  in  yaw 

B.  =  Coefficients  of  hydrodynamic  forces  in  sway 

C.  =  Coefficients  of  hydrodynamic  forces  in  surge 
i  =  1,2,3,4 

N.,  N~,  Y  ,  X.  are  coefficients  of  aerodynamic  forces 

and  moments 

Aerodynamic  forces  and  moment 

xa 

X    '    =  = =  X.  *    U    '    =   X      ■    cos   \ba 

A         |^L  £2   Ua2  i       a  ai 


N  Na' 

N.  ■    =  =  N'u'v'+N'v'    =  -x-  sin   2i|;a+Na;    sin  ipa 

A  |a   £3    Ua2  2  2 


Y 

-  Y    •   V    ■    -  Y      '    sin  tya 


A  I3.  £2   Ua2  la  &1 
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Which 


X    '    -  Xa    '    -   -0.015 
1  ax 


y    i    =    y       f    =    -0.056 
1  ax 


N, '    =   N       '    =    -0.0046,         N    '    =   -N       '    =   -0.017 
1  al  2  a2 
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Figure  5.   Variation  of  Yaw  Versus  Time  for  Wind  Direction 
000° ,  Velocity  1  mile/hr. 
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Figure  6 


Variation  of  Sway  Versus  Time  for  Wind  Direction 
000°,  Velocity  1  mile/hr. 
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Figure  7.   Variation  of  Yaw  Versus  Time  for  Wind  Direction 
000° ,  Velocity  3  miles/hr.   ' 


60 


•o 


a. 
o 


2> 

rs 

\ 

*- 

< 

3 

•>■» 

CO 

'O 

*> 

I 

I 


1 

)0O 

loo 

lOO 

too 

CfOC 

1 

I 

i 
i 

v_ 

T/Alt  (Sec) 


Figure  8.   Variation  of  Sway  Versus  Time  for  Wind  Direction 
'    000°,  Velocity  3  miles/hr. 


61 


/ 


% 

ICO 

ic» 

gft 

^i^Do 

d£'pd 

* 

I 

1 

a.  » 
> 

/^ 

«0 

I 

if 

\ 

1 

~ft/JE  C  S«<.) 


Figure  9.   Variation  of  Yaw  Versus  Time  for  Wind  Direction 


000 


Velocity  5  miles/hr. 
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Figure  10. 


Variation  of  Sway  Versus  Time  for  Wind  Direction 
000° ,  Velocity  5  miles/hr. 
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Figure  11.   Variation  of  Yaw  Versus  Time  for  Wind  Direction 
090° ,  Velocity  3  miles/hr. 
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Figure  12.   Variation  of  Sway  Versus  Time  for  Wind  Direction 
090° ,  Velocity  3  miles/hr.  *  ' 
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Figure  13.   Variation  of  Yaw  Versus  Time  for  Wind  Direction 
090°,  Velocity  5  miles/hr. 
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Figure  14.   Variation  of  Sway  Versus  Time  for  Wind  Direction 
090°,  Velocity  5  miles/hr. 
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Figure  15.   Variation  of  Yaw  Versus  Time  for  Wind  Direction 
090°,  Velocity  7  miles/hr. 
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Figure  16.   Variation  of  Sway  Versus  Time  for  Wind  Direction 
"  090°,  Velocity  7  miles/hr. 
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Figure    17 


lSO^v^V"?  V6fS^  Time  for  Wind  Direction 
-lou  ,  Velocity  1  mile/hr. 
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Figure  18.   Variation  of  Sway  Versus  Time  for  Wind  Direction 
180°,  Velocity  1  mile/hr. 
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Figure  19 


Y?^ation  of  Yaw  Versus  Time  for  Wind  Direction 
180° ,  Velocity  3  miles/hr. 
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Figure  20.   Variation  of  Sway  Versus  Time  for  Wind  Direction 
180°,  Velocity  3  miles/hr. 
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Figure    21 


Y«^ation  °f  Yaw  VersVis  Time  for  Wind  Direction 
180°,  Velocity  5  miles/hr. 
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Figure  22.   Variation  of  Sway  Versus  Time  for  Wind  Direction 
180° ,  Velocity  5  miles/hr. 
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Figure  23. 


Variation  of  Yaw  Versus  Time  for  Wind  Direction 
180° ,  Velocity  7  miles/hr. 
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Figure  24. 


Variation  of  Sway  Versus'  Time  for  Wind  Direction 
180°,  Velocity  7  miles/hr. 
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Figure  25. 


Variation  of  Yaw  Versus  Time  for  Wind  Direction 
090°,  Velocity  5  miles/hr,  Rudder  Deflection  15° 
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Figure  26 . 


Variation  of  Sway  Versus  Time  for  Wind  Direction 
090° ,  Velocity  5  miles/hr,  Rudder  Deflection  15' 
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Figure  27.   Variation  of  Yaw  Versus  Time  for  Wind  Direction 
180°,  Velocity  3  miles/hr,  Rudder  Deflection  15 
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Figure  28. 


Variation  of  Sway  Versus  Time  for  Wind  Direction 
180° ,  Velocity  3  miles/hr,  Rudder  Deflection  15° 
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Figure   29. 


Variation  of  Yaw  Versus  Time  for  Wind  Direction 
180° ,  Velocity  5  miles/hr,  Rudder  Deflection  15° 
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Figure  30. 


Variation  of  Sway  Versus  Time  for  Wind  Direction 
180°,  Velocity  5  miles/hr,  Rudder  Deflection  15° 
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Figure  31.   Variation  of  Yaw  Versus  Time  for  Wind  Direction 
000° ,  Velocity  1  mile/hr,  Automatically  Steered. 
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Figure  32.   Variation  of  Sway  Versus  Time  for  Wind  Direction 
000°,  Velocity  1  mile/hr,  Automatically  Steered. 
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Figure  33.   Variation  of  Yaw  Versus  Time  for  Wind  Direction 
000°,  Velocity  3  miles/hr,  Automatically  Steered 
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Figure  34. 


Variation  of  Sway  Versus  Time  for  Wind  Direction 
000°,  Velocity  3  miles/hr,  Automatically  Steered 
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Figure  35.   Variation  of  Yaw  Versus  Time  for  Wind  Direction 
000° ,  Velocity  5  miles/hr,  Automatically  Steered 
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Figure  36.   Variation  of  Yaw  Versus  Time  for  Wind  Direction 
090°,  Velocity  3  miles/hr,  Automatically  Steered 
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Figure  37.   Variation  of  Sway  Versus  Time  for  Wind  Direction 
090°,  Velocity  3  miles/hr,  Automatically  Steered. 
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Figure  38.   Variation  of  Yaw  Versus  Time  for  Wind  Direction 


180 


Velocity  1  mile/hr,  Automatically  Steered. 
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Figure  39.   Variation  of  Sway  Versus  Time  for  Wind  Direction 
180° ,  Velocity  1  mile/hr,  Automatically  Steered. 
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Figure  40. 


Yaw  Gain  Constant  Versus  Yaw  Rate  Gain  Constant 
for  Wind  Direction  000° ,  Velocity  1  mile/hr. 
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Figure  41. 


Yaw  Gain  Constant  Versus  Yaw  Rate  Gain  Constant 
for  Wind  Direction  000° ,  Velocity  3  miles/hr. 
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Figure  42. 


Yaw  Gain  Constant  Versus  Yaw  Rate  Gain  Constant 
for  Wind  Direction  000° ,  Velocity  5  miles/hr. 
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Figure  43.   Yaw  Gain  Constant  Versus  Yaw  Rate  Gain  Constant 
for  Wind  Direction  090°,  Velocity  3  miles/hr. 
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Figure  44. 


Yaw  Gain  Constant  Versus  Yaw  Rate  Gain  Constant 
for  Wind  Direction  180° ,  Velocity  1  mile/hr. 
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Figure    45 


lA/M/D    VELOC'li    (milt) 
Real   o^  Versus   Wind  Velocity   for  Head  Wind, 
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Figure    46.      Real   a.    Versus  Wind  Velocity   for   Beam  Wind 
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Figure  47.   Real  o.  Versus  Wind  Velocity  for  Stern  Wind 
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Figure  48. 


Real  a   Versus  Wind  Direction  for  7  miles/hr  Wind 
Velocity. 
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Figure  49 


Real  a   Versus  Wind  Direction  for  9  miles/hr  Wind 
Velocity. 
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